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Standing Electromagnetic Solitons in Degenerate Relativistic Plasmas 
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The existence of standing high frequency electromagnetic (EM) solitons in a fully degenerate 
overdense electron plasma is studied applying relativistic hydrodynamics and Maxwell equations. 
The stable soliton solutions are found in both relativistic and nonrelativistic degenerate plasmas. 


A significant amount of recent publications describe electromagnetic (EM) waves in relativistic plasmas and majority 
of them discuss possible roles of these waves in different astrophysical phenomena. Highly relativistic plasmas are 
observed in the cores of white dwarfs Q, in magnetosphere of pulsars [l|, in the MeV epoch of the early Universe 
fH and additionally, they probably show up in the bipolar jets in Active Galactic Nuclei (AGN) [3]. Plasma can 
be relativistic in two following cases: either bulk velocities of fluid cells should be close to the speed of light, or the 
kinetic energy of particles should be greater then their rest energy. In compact objects, such as white dwarfs and 
magnetars, the number densities of electrons is believed to be roughly between 10 26 cm _3 and 10 34 cm -3 Q. High 
density plasma can be produced in the laboratory as well, indeed contemporary petawatt laser systems have the focal 
intensities I = 2 x 10 22 W/cm 2 [71]. Moreover, pulses with higher than I = 10 26 W/cm 2 intensities are expected to be 
achieved soon Q. Superdense plasmas might be formed with densities in the range of 10 23 cm -3 and 10 28 cm -3 [§], 
during the interaction of such EM pulses with solid or gaseous targets. Such plasma will be opaque for conventional 
laser systems operate at wavelengths A ~ lpm. The Linac Coherent Light Source (LCLS) is an X-ray free-electron 
laser produce femtosecond powerful pulses of coherent soft and hard X-rays with wavelengths from 2.2 nm to 0.06nm 
0 - Exploiting the possibility to focus X-ray laser beams on a spot with down to laser wavelength, the focal intensities 
I ~ 7 x 10 25 W/cm 2 are expected to be reached [ll| . Successful operation of X-ray free-electron lasers in different 
centers world wide E3 opens up new perspectives to study the EM pulse penetration and its subsequent dynamics 
in super dense plasma in laboratory conditions. 

Highly compressed plasma with an average interparticle distance smaller than their thermal de Broglie wavelength, 
can be considered as a degenerate Fermi gas. When plasma density increases, the more ideal it becomes and the 
interactions of its particles can be neglected [l3]. 

EM solitons in classical relativistic plasma is being studied intensively Q but existence and stability of solitary 
solutions in degenerate quantum relativistic plasma are investigated mostly for low frequencies (see fl5j and references 
therein). The publication goal is to consider existence of a standing, high frequency EM soliton in the relativistic 
degenerate electron plasma. Importance of the standing soliton solutions for overall dynamics of EM pulses is estab¬ 


lished theoretically [If 


El as well as experimentally for classical relativistic plasma [20|. These publications state, 
that during interaction of a circularly polarized strong laser pulse and a plasma, part of the laser energy is trapped 
in non-propagating soliton-like pulses. Similar dynamics is expected in the case of strong EM pulse interaction with 
degenerate electron plasma. 

Plasma can be considered cold, if the thermal energy of its electrons is negligible compared to their Fermi energy. 
In this case temperature can be assumed to be zero, even though it is of the order of 10 9 X [2l| . For the Fermi energy 

of electrons we have €f = m e c 2 (l + I? 2 ) — 1 , where R = pp/m e c and pf is the Fermi momentum. The latter is 

related to the proper density of electrons n by the following equation pp = m e c (n/ric ) 1 ^ 3 , where n c = 5.9 x 10 29 cm -3 
is the normalizing critical number-density |22[. Therefore, when n > n c , electrons move with relativistic momentum 
inside plasma unit cells and the plasma can be told as relativistic. 

Our investigation is based on the Maxwell equations and fluid model of relativistic electron plasma. The ions are 
considered to form a stationary neutralizing background. We begin with the manifestly covariant form of the fluid 
equations for electrons 


QjyaP 

dxP 


eF^nUp 


(1) 


here d a = d/dx a = (c 1 d/dt, V); the Greek indexes take values from 0 to 3; T a P is the energy-momentum tensor 
describing the plasma electrons with charge —e, mass m e and the proper density n; the metric tensor is g a @ = 
diag ( 1,-1,—1,-1); U a = ( 7,7 V/c) denotes the local four velocity, here 7 = (l — V 2 /c 2 ) 1 ^ 2 ; ( U a U a = 1). This 
equation implies the conservation of energy and momentum. The change of momentum through the collisions is 
neglected. 

We assume, that the total number of electrons is conserved, thus the following continuity equation is held 





2 


dnU a 

dx a 


( 2 ) 


The EM field can be expressed through a tensor F a @ = [E,B]. The Maxwell equations in these notations are 
dpF a P = —(47t/c)J q and e a ^ s dpF 1 s = 0. Here J a = (cp, J), where J is the current density and p is the total 
charge density of the plasma. 

We use the energy momentum tensor of ideal isotropic fluid T a P = wU a U& — g a/3 P, where w = E + P is the 
enthalpy per unit volume, P is the pressure and E is density of the rest frame internal energy. If nT/P « 1, plasma 
can be treated as completely degenerate Fermi gas and the following equations are satisfied [23] 


P = 


m^c 5 
3ir 2 H 3 


f(R) 


( 3 ) 


E = 


mfc 5 

3n 2 h 3 


R 3 (1 + R 2 ) 1/2 


f(R) 


( 4 ) 


where 


8 / {R) = 3 sinlV 1 R + R.( 1 + R 2 ) 1/2 (2 R 2 - 3) (5) 

The equation of state for the degenerate gas is P oc n r , with T = 5/3 for nonrelativistic case (R « 1) and T = 4/3 
for ultrarelativistic case (R » 1). 

The model of plasma described above implies that the electron distribution function remains locally Juttner- 
Fermian. In case of zero temperature this results in thermodynamical quantities, depending only on density E(n), P(n) 
and w(n). Of course, these quantities are functions of x a through the relation n = _/V/ 7 , here N is the electron density 
in the laboratory reference frame. The considered system is isentropic (furthermore, as temperature approaches zero, 
entropy tends to zero too). Hence, the following thermodynamic equality is held d(w/n) = dP/n and taking into 
account this thermodynamic equality and making some standard manipulations (e.g. |2J|), Eq. © can be represented 
in the form of the following system 


- (Gp) + ?n e c 2 V (G 7 ) = —eE+ [V x n] 


( 6 ) 


—ft = vx [v x n\ 

dt 


( 7 ) 


here, for the generalized vorticity we have = — (e/c) B + Vx (Gp) , where p = 7 m e V denotes the hydrodynamic 
momentum; G = G(n) can be called the density dependent ’’effective mass” factor of electrons G = w/m e nc 2 = 

1/0 

(l + R 2 ) . Now dynamics of the degenerate plasma can be completely described by Eqs. ©-© together with 

Continuity and Maxwell equations. In other words, the set of equations is complete. The analogous set of equations 
is derived in Ref. [24] for classical relativistic plasma obeying Maxwell-Juttner statistics, where G is a function of 

temperature G = G(T). In contrast, for degenerate plasma w/m e nc 2 = (l + R 2 ) " and as a result the effective mass 


1 + (n/n c ) 2/3 


11/2 


holds for 


factor of electrons depends just on their proper density. The corresponding relation G = 

any ratio n/n c , thus for any strength of relativity [25j -[26]. 

We make use of the expressions for fields B = V x A and E = —(l/c)dA/dt — 'S/ip where A and ip are vector and 
scalar potentials respectively. Applying the Coulomb gauge condition V ■ A = 0, the Maxwell equations take the 
following form: 


d 2 A 

~8P 


d 

AA+c— (V<p) — 47 tcJ = 0 


( 8 ) 


and 
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A ip = —4:TTp (9) 

where J = —eynV, is the current density and p = e(no — jn) is the charge density, with no denoting electron (ion) 
equilibrium density. We use equations © and © to describe wave dynamics in unmagnetized plasma. Eq. © makes 
clear that if generalized vorticity fi was zero everywhere once, it will stay zero always. Therefore Eq. © will reduce 
to 


(Gp-^A) + V (m e c 2 G'y- e<p) = 0 (10) 

Our goal is to find one dimensional localized solutions for equations Let us assume, that every variable 

depends on nothing but coordinate z and time t. As transverse component of gradient is zero, Eq. (1101) easily gives 
px = eAj_/(cG). Integration constant is zero, because p should be zero at infinity, where fields vanish. Coulomb 
gauge condition requires A z = 0, thus the longitudinal motion of the plasma is driven just by the ” ponderomotive” 

*1 / o 

pressure (~ pj_) acting via the relativistic 7 factor in Eq. ©3 (7 = [1 + (pi +pl) /m 2 c 2 ] ). The EM pressure 

forces electrons to move in z direction, the plasma density changes and charge separation occurs. Longitudinal motion 
of the plasma is described by the following set of equations: 


o o 

Gp z + Try (m e c 2 G 7 - etp) = 0 


dt ~ dz 

while the continuity © and Poisson’s equations © become 


(ID 


| 7 „+^(n 7 U) = 0 


( 12 ) 


d 2 <p 

IhT 1 


= 47re(n7 — no) 


(13) 


The transverse component of the current density is J^= (?re 2 /cG) Aj_ and substituting it into Eq. ©, we get 


d 2 A 


dt 2 


± 2 

— c 


d 2 A± | n2 / n_Go 
3z 2 e \n 0 G 


A_l = 0 


1 j 2 

where no is electron (ion) equilibrium density and f l e = (47re 2 no/nr*) is the Langmiur frequency of the electron 

1 /2 

plasma. In this expression m* denotes effective mass of electron m* = m e Go , where Go = [l + Rq] and Rq = 
{no/n c ) 1/3 . 

We are searching for solitary stationary solutions of Eqs. (111I) - (I14[) . Assuming electromagnetic wave is circularly 
polarized, the vector potential can be rewritten as follows Ax = [A{z) cos(wt), A(z) sin(wt)]. Here the amplitude 
A{z) is a real valued function depending only on coordinate z and nothing else, w denotes the frequency. It is now 
convenient to introduce the following dimensionless quantities: T = e<p/ ( m e c 2 Go ), a = eA/ (c 2 m e Go), n = n/no, 
t = fl e t, z = (f2 e /c) 2 . In our case p z = 0 and integrating Eq. (fill) we get the relation 


(14) 


Ty = 'I' + 1 


(15) 


where 


r = — = 

G 0 


A + RqU 2 / 3 
1 + R o 


1/2 

and 7 = (l + a 2 /T 2 ) . Using Eqs. (115l) - (fl6l) we obtain the following relationships: 


(16) 
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n = 


Rl L 


(1 + ^) (1 + ^) 2 — a 2 


- 1 


3/2 


(17) 


7 = 


(1 + 4/) 


^/(l + tf) 2 -a 2 

while for the electron density in laboratory frame A/" (= ny) we have 


(18) 


AT 


(1+Vh) 


, 3/2 


(1 + 4>) 2 - a 2 


1/2 


(1 + 40 2 


(1 +a 2 



(19) 


where e = Rq/ (l + i?p). 

The Eqs. (fl3l) and (1141) now can be reduced to the following set of ordinary differential equations: 


d 2 a 

te 2 


(1 — w 2 )a + 


AT (a, dq e) \ 

1 + T ) 


a = 0 


( 20 ) 


+ 1 — AA (a, 4/, e) = 0 (21) 

The coupled system of nonlinear equations (fTPl) - (ETl) . where AT (a, 4/, e) is defined by Eq. (fTPl) describes the structure 
of circularly polarized transverse and longitudinal localized fields in degenerate plasma for arbitrary values of Rq 
(0 < e < 1). The general solution of Eqs. (fTPl) - (l2Tl) cannot be expressed in terms of elementary functions, it requires 
numerical methods. However, these equations can be integrated analytically for certain limiting cases. Taking into 
account that for \z\ -+ oo the fields vanish a —> 0, 4/ —> 0 and AT —> 1 we conclude that the standing solitary solutions 
may exist in an overdense plasma when the frequency of the soliton is less than the electron plasma frequency, to < 1 
(w < - in dimensional units). For any fixed plasma density (i.e. fixed e) the frequency u> is the only parameter 

in Eqs. m-m that defines the characteristics of the solitary solutions. For slightly overdense plasma u> < 1 the 
fields are small a <C 1, 4 1 <C 1 while the characteristic width of the structure L~ (1 — w 2 ) -1 / 2 >> 1. Applying the 
quasi-neutrality condition of the plasma (AT «1) from Eqs. (fTHl) and (1211) we get 4/ ~ (n 2 /2) (1 — e/3) and now Eq. 
(EDI) reduces to the simple equation 


d 2 a 

dz 2 


(1 - ui 2 )a + - (1 


e/3) a 3 = 0 


( 22 ) 


with the well know soliton solution 


CL — 2 a 


' 1-w 2 
1 — e/3 


Sech 


:\/l — uj 2 


(23) 


Analytical solutions of Eqs. m-m can be obtained even for an arbitrary amplitude fields provided that plasma 
density is small no/n c « 1 (Rq « 1), implying that the degenerate electron gas is not relativistic. At this 
end we would like to emphasize that our consideration is valid if the average kinetic energy of electrons (~ €f) 
is larger than their interaction energy (~ e 2 n^ 3 ). This condition is fulfilled for a sufficiently dense plasma when 

no >> ^2 m e e 2 / (37r 2 ) 2 ^ 3 h 2 ^ = 6.3 x 10 22 cm -3 . For the densities no in the range from 10 24 cm -3 to 10 28 cm -3 the 

corresponding parameter e is in the range from 1.4 x 10 -4 to 6.2 x 10~ 2 . Since e is rather small for the nonrelativistic 
densities, we can safely assume that e —> 0 in Eqs. m-m and then obtain the following relations: 
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FIG. 1: The profiles of EM field potentials a, VP and the electron density N in case the ultrarelativistic plasma e = 1 for the 
different soliton frequencies (a) to = 0.8 and (b) uj = 0.7683. 


= 1 - \fl 


(24) 


the wave equation (12T1) now reduces to 


N=1 


dz 2 


\/l + . 


(25) 


dz 2 


-\~ UJ CL — 




1 + 


dz 2 


vT 


= 0 


(26) 


Similar to (El- ES system equations have been derived solved in Refs. [27|, [28| where existence and stability of 
standing EM pulses in cold but classical electron plasma were studied. The soliton solution of Eq. (l26l) is found to be 


2\/l — w 2 cosh (zjl — uj 2 ) 

° = cosh 2 (zVl^)+u, 2 -l (27) 

The amplitude of this single hump soliton a m (= a(0)) depends on the soliton frequency by the relation a m = 
2 VT — u) 2 /w 4 . The electrons are depleted from the region of pulse localization. The minimal density at the center 
of soliton is J\f(z = 0) = 1 — 4 (l — w 2 ) /w 4 [Ii|. Note that for u 4 1 the amplitude of soliton decreases and the 
solution (1271) coincides with Eq. C3l) where e « 1. With a decrease of w the soliton amplitude increases while the 
corresponding density well deepens. For u> = u> c = \f2/2> = 0.816 the electron cavitation takes place Af(z = 0) = 0 
while the amplitude of the soliton attains its maximal value a m = JS. For a m > V?> the solution contains a region 
where the electron density is negative which implies that wave breaking takes place. Thus, in the nonrelativistic 
degenerate electron plasma the standing solitons exist provided ui c < to < 1. Similar conclusion has been made 
in Ref. f29l| by numerical analysis of the coupled set of nonlinear Schrodinger-Poisson equations, where the Bohrn 
potential was introduced in the system. In the fluid approach the Bohm potential is related to the ’’quantum force” 
due to electron tunneling. Our estimations show that this force is small and even for w = uj c it does not change the 
solution EH qualitatively. 

Now we consider the case of ultrarelativistic degenerate plasma Rq » 1 (e —;> 1). The system (THHl - dTIl) reduces to 
the following set of equations 


d 2 a 

dz 2 


+ co 2 a + a (a 2 — (1 + 4/) 2 ^ = 0 


(28) 


and 


^ + 1 + (1 + T) (a 2 - (1 + 4-) 2 ) = 0 


(29) 
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FIG. 2: The profiles of EM field potentials and the electron density in case the relativistic plasma with e = 0.5 for the different 
soliton frequencies (a) u = 0.8 and (b) lj = 0.7694. 


This system of equations has a first integral (a Hamiltonian) 


'(N>' 

2 

da 

dz 


dz 


- uj 2 a 2 + 2T - 



(l + 'F ) 1 2 



(30) 


where the zero boundary conditions for the fields 'F, a and their derivatives at \z\ —> oo have been used. 

Fig.l shows the numerical solution of the system for different values of the soliton frequency w. One can see 

that the structure of the solitary solutions in ultrarelativistic plasma is similar to one, obtained in the nonrelativistic 
case. The solution has a continuous spectrum with co c < ui < 1 where w c = 0.7683, and it is composed of a single 
maximum of field potentials (a, H/) and depleted electron density (AT) region at the center of structure. For w —> 1 the 
numerical solution coincides with Eq. (1231) (using e = 1). For lo ~ w c , amplitude of the soliton becomes relativistically 
strong a m ~ 2.267, while at the center of the soliton Af ~ 0, i.e. the electron cavitation takes place. 

Numerical analysis demonstrates that equations (l20l) - (l2lT) have soliton-like solutions for any finite value of Rq. As 
an example in Fig.2 we plot the numerical solutions of the system for R 0 = 1 (e = 0.5). The soliton exists for 
w c < w < 1, where w c = 0.7694 and at this frequency a m — 2.162. 

To verify the stability of obtained solutions, we performed ID simulations applying the numerical code developed in 
Ref. [01 . The simulations demonstrate exceptional stability of the solution for entire rage of soliton existence described 
above. 

To summarize, in this paper we considered circularly polarized high frequency EM solitons in degenerate electron 
plasma. We used Maxwell and relativistic fluid equations to demonstrate possibility of stable existence of solitons in 
overdense plasma (w < fi e ). Soliton exists for entire range of physically allowed electron densities, presumably for 
no = 10 24 cto ~ 3 4 5 6 7 8 9 10 and higher, i.e. in both relativistic and nonrelativistic degenerate plasmas. Intensity of the solitons 
can be small for w —> fi e and becomes relativistically strong (a m > 1) for u> —> lu c . Appears to be, that cavitation 
of plasma takes place in both, relativistic and nonrelativistic degenerate plasmas. The model described above can be 
straightforwardly generalized for underdense plasma. This is not done in the scope of this paper deliberately. 

The present results could have rather interesting implications to describe and understand A'-ray pulses emerging 
from compact astrophysical objects. It is also of particular interest to describe interaction of intense laser pulses and 
dense degenerate plasma, as such experiments are becoming feasible with the new generation of intense lasers. 
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